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ABSTRACT 

A model for irradiated dust disks around Herbig Ae stars is proposed. The model is based on the 
flaring disk model of Chiang & Goldreich (1997, henceforth CG97), but with the central regions of the 
disk removed. The inner rim of the disk is puffed up and is much hotter than the rest of the disk, because 
it is directly exposed to the stellar flux. If located at the dust evaporation radius, its reemittcd flux 
produces a conspicuous bump in the SED which peaks at 2-3 micron. We propose that this emission 
is the explanation for the near-infrared bump observed in the SEDs of Herbig Ae stars. We study for 
which stellar parameters this bump would be observable, and find that it is the case for Herbig Ae stellar 
parameters but not for T-Tauri stars, confirming what is found from the observations. We also study 
the effects of the shadow cast by the inner rim over the rest of the flaring disk. The shadowed region 
can be quite large, and under some circumstances the entire disk may lie in the shadow. This shadowed 
region will be much cooler than an unshadowed flaring disk, since its only heating sources are radial 
radiative diffusion and possible indirect sources of irradiation. Under certain special circumstances the 
shadowing effect can suppress, or even completely eliminate, the 10 micron emission feature from the 
spectrum, which might explain the anomalous SEDs of some isolated Herbig Ae stars in the sample of 
Meeus et al. (2001). At much larger radii the disk emerges from the shadow, and continues as a flaring 
disk towards the outer edge. The emission from the inner rim contributes significantly to the irradiation 
of this flaring disk. The complete semi-analytical model, including structure of the inner edge, shadowed 
region and the flared outer part, is described in detail in this paper, and we show examples of the general 
behavior of the model for varying parameters. 



1. INTRODUCTION 

The hypothesis that Herbig Ae/Be stars (HAeBe; Her- 
big 1960) are intermediate mass pre-main-sequence stars 
is today generally accepted. The large infrared (IR) 
excess observed from these stars can be naturally ex- 
plained in this context by emission from a protostel- 
lar/protoplanetary disk. This picture has been given more 
credibility in recent years by interferometric observations 
at millimeter wavelengths (Mannings & Sargent 1997, 
2000). Over the years, various disk models have been pro- 
posed for pre-main-sequence stars, among which are disks 
powered by accretion (Lin & Papaloizou 1980, Bell & Lin 
1994), irradiated non-accreting (passive) disks (Kenyon 
& Hartmann 1987, Chiang & Goldreich 1997 henceforth 
CG97), and irradiated accretion disks (Ruden & Pollack 
1991, Calvet et al. 1991, 1992, D'Alessio et al. 1998, 1999, 
Bell 1999). The spectral energy distributions (SEDs) of 
Herbig Ae stars seem to be best fitted by passive irradiated 
flaring disks, at least at mid- and far-infrared wavelengths 
(Chiang et al. 2001). At near-IR wavelengths, however, 



the flaring disk picture fails to agree with the observed 
SEDs. The same problem is encountered when fitting non- 
irradiated accretion disk models to the SEDs. Hillenbrand 
et al. (1992) noticed that in 29 out of 51 HAeBe stars (their 
Group I sources) the SED at near and mid-IR wavelengths 
was well fit by models of geometrically flat, optically thick 
disks with very large accretion rates, but that these disks 
needed to be optically thin in their inner parts, to account 
for the distinctive inflection shown by all HAeBe stars at 
1.2 < A < 2.2 fxm (see also Lada & Adams 1992). Bet- 
ter measured SEDs for a number of Herbig AeBe stars of 
spectral type A (HAe in the following) have confirmed the 
Hillenbrand et al. (1992) results. 

All the observed isolated HAe stars have large near- 
infrared emission in excess of the photospheric one with 
very similar dependence on wavelength, that can be de- 
scribed as a near-infrared "bump" : the excess monochro- 
matic luminosity v¥ v is negligible around 1 /im, peaks 
around 2 /im and decreases slowly with wavelength to 
about 8/xm, where the silicate emission feature starts 
(Meeus et al. 2001; Natta et al. 2001). These stars 
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have very little extinction, so that reddening in the near- 
infrared is negligible. 

The negligible opacity of the inner disk region, required 
to explain the near-infrared bump is not consistent with 
the high accretion rates necessary to account for the large 
near-infrared excess, as discussed by Hartmann, Kenyon 
& Calvet (1993). These authors suggested that, given 
the difficulties of applying disk models to HAeBe stars, 
one should consider the hypothesis that the emission orig- 
inates in a dusty circumstellar envelope, possibly contain- 
ing very small grains transiently heated by ultraviolet pho- 
tons. This possibility, however, has been ruled out by 
theoretical models (Natta & Krugel 1995) and by ISO ob- 
servations (Meeus et al. 2001). 

An alternative explanation for the near-infrared be- 
haviour of HAe stars has recently been proposed by Natta 
et al. (2001). Their model starts from the notion that 
the mid- and far-infrared excess of HAe stars can be well 
fitted with a passive irradiated flaring disk, as was men- 
tioned above. As expected, these models fail to explain 
the observations in the near-infrared. However, Natta et 
al. have noticed that, if one truncates the flaring disk at 
a certain radius (i.e. introduce a hole), the emission from 
the inner rim of this disk is non-negligible. This inner 
rim will have a considerable covering fraction around the 
star, and it will acquire a temperature much hotter than 
what the flaring disk interior at the same radius would be. 
This is because the flaring disk receives the stellar radi- 
ation under a small grazing angle, and will therefore be 
relatively cool, while the inner rim is in full sight of the 
star and will therefore be hot. In fact, if one assumes that 
the inner edge of the disk is naturally located at the radius 
where dust grains acquire their evaporation temperature, 
then the emission from this inner rim is going to produce 
a near-infrared bump very reminiscent of the observed one 
in shape and intensity. 

In this paper we wish to proceed further along this path, 
and investigate a number of issues that have not been ad- 
dressed by Natta et al. (2001). Firstly, we will formalize 
the description of the structure of a disk with a large inner 
radius, using the flaring disk models of CG97 as a start- 
ing point. In particular, we wish to investigate how the 
presence of the hot inner rim modifies structure and ap- 
pearance of the rest of the disk. Since the inner rim is 
much hotter than the CG97 disk at the same radius, its 
vertical scale height will be larger as well. In other words, 
the inner rim of the disk will puff-up. This will enhance 
the near-infrared emission, but also cast a shadow over the 
disk behind it. Part of the disk which would normally be 
directly illuminated by the stellar radiation, will now be 
in the shadow, and therefore collapse and cease to emit 
IR radiation. At large enough radii, it is possible that 
the flaring shape of the CG97 disks will be conserved, and 
that the disk will grow out of the shadow, unchanged by 
the presence of the rim. We thus propose a scenario in 
which the emission from the reprocessing disk comes from 
a hot inner rim, caused by dust evaporation, plus the non- 
shadowed part of a truncated CG97 disk. 

A second point we will discuss concerns the behaviour of 
the near-infrared emission as function of the spectral type 
of the star. The IR spectra of T Tauri stars (TTS) are well 
described by disk models with small or zero inner holes 



(Beckwith et al. 1990; Kenyon, Yi & Hartmann 1996) and 
moderate accretion rates (Gullbring et al. 1998). The few 
known prc-main-scquencc stars of spectral type F and G 
seem closer in properties to TTS than to HAe stars (Hillen- 
brand et al. 1992; Hartmann et al. 1993). We show that 
this behaviour can be explained by the same disk models 
and that it is caused by a combination of effects. First 
of all, the relative importance of the NIR bump emission 
from the inner edge in comparison to the disk emission 
is smaller for stars of later spectral type. Secondly, the 
stellar atmospheric spectrum shifts towards longer wave- 
lengths as the effective temperature of the star decreases, 
swamping the rim emission. Finally, it is possible that 
TTS have on average higher accretion rates, and therefore 
more optically thick inner disks, than HAe stars. 

The structure of the paper is as follows. Section 2 dis- 
cusses the equations, Section 3 describes how the SED is 
computed. The results will be presented in Section 4, and 
we conclude with a discussion in Section 5. 

2. THE STRUCTURE OF A DISK WITH A LARGE INNER 

HOLE 

In this section we will derive the equations governing 
the structure of a passively irradiated circumstellar disk 
with an inner hole. The main effect of the inner hole on 
the disk structure is a vertical boundary which is irradi- 
ated directly by the star, heating it to temperatures higher 
than what would be expected at the same location for a 
thin disk with grazing irradiation only. Because of the 
higher temperature, the disk boundary will be puffed up, 
and the material immediately outside the inner rim will be 
shadowed, leading to a decreasing scale height. At large 
distances from the star it is possible that the disk becomes 
high enough to emerge from the shadow of the inner rim, 
in which case it may continue to flare like a disk without 
a puffed-up inner rim. However, it is not certain that the 
disk will be able to leave the shadow. This is a dynami- 
cal question and the answer may depend upon the history 
of the disk and many other parameters (e.g. was the disk 
flaring when the inner boundary started to emerge and to 
puff up?). 

Figure Bl shows our scenario, where we assume that the 
flared part of the disk exists when it can. 

We will divide the disk structure calculations in three 
parts (inner rim, shadowed region, and flared component) 
and derive the basic equations governing the structure in 
the three regions. We start with the well-known flared 
part. 

2.1. Flaring disk 

For the the flaring part of the disk we follow the line 
of Chiang & Goldreich (CG97) with some minor modifi- 
cations. For completeness, and in order to introduce the 
nomenclature we repeat the basic equations without phys- 
ical discussion. 

2.1.1. Disk interior 

Direct stellar radiation impinges onto the disk at an an- 
gle a given by a = QAR*/R + Rd(H cg /R)/dR, where R* 
is the radius of the star, R the distance from the star in 
the disk plane and H cg the height of the disk photosphere 
above the midplane. Anticipating a nearly powerlaw be- 
havior, we replace Rd/dR by 7 — 1, where 7 is a slowly 



3 



varying function of the order of 9/7. We can thus rewrite 
the expression for a as: 



0.4R, i\^c g 



(1) 



The (7 — 1) factor is self-consistently determined as a func- 
tion of R, using a numerical method devised by Chiang et 
al. (2001). In this way energy conservation is ensured, 
while the assumption of a constant (7 — 1) can cause de- 
viations by a factor <~ 2. The stellar flux impinging with 
this angle into the disk is: 
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This equation assumes that the entire star is visible from 
the considered point on the disk surface. If the disk would 
stretch all the way to the star, the bottom half of the star 
would not be visible and a correction would have to be 
made. However, since we consider disks with large inner 
holes, the transition from the star being invisible to the 
star being fully visible will be sharp, and we ignore this 
effect. 

The flux F ; * r is absorbed in the uppers layer of the disk, 
which will re-radiate half of the flux away from the disk, 
and half down into the disk's deeper layers. A fraction 
tp s of this flux will be absorbed by the interior, where ip s 
is a dimensionless number in between and 1 that ac- 
counts for the possibility that the disk interior is not fully 
optically thick to the emission of the surface layer. We 
assume the disk interior to be isothermal with a tempera- 
ture T;, which is a very good approximation for optically 
thick disks. This interior will emit a flux F em n given by 
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where ipi accounts for the possibility that the disk interior 
is not fully optically thick to its own emission. For a fully 
optically thick disk one has ip s = ipi = 1. The formulas for 
determining ip s and ipi for not optically thick cases will be 
given in appendix A.l. 

By equating (l/2)^ s F ; * r = -F C mit, one can solve for T; 
and find: 
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The disk will be almost isothermal in the vertical di- 
rection, leading to a Gaussian density distribution with a 
constant pressure scale height h cg given by 
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with T c = GM*/j,m p /kRi, is the virial temperature at the 
stellar surface. The ratio of the disk surface height H cg 
to the pressure scale height h cg is a dimensionless number 
Xcg of order unity: 
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Usually this Xcg will lie somewhere in between 2 and 6. 
But it can be self-consistently determined from the graz- 
ing angle a, the surface density S and the Planck mean 
opacity at the stellar temperature Kp(T^). The way to do 
this will be described in appendix A. 2. 



2.1.2. Disk surface layer 

Since the disk surface layer is defined as the layer of 
matter in direct sight of the central star up to grazing op- 
tical depth r a ~ 1, the temperature in this layer can be 
estimated from the optically thin expression: 

where e s is the ratio of the Planck mean opacities at T s 
and T». Quantities describing the surface layer carry the 
subscript s. The Planck mean opacity is defined as 
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The flux from the surface layer in both the upwards and 
the downwards direction is: 
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F s = 2ttAE / B v (T s )n v dv 
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= 2AZk p (T s )<jT s 4 , 

where AE is the surface density of the layer. Energy con- 
servation implies F s = (ipi /ip s )<rT^, so the surface density 
of the layer is: 

1 iPi 
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2k p (T s ) Vs T s 4 ' 
In our treatment of the surface layer we ignore scatter- 
ing. It should be pointed out that, while the effect of this 
approximation on the structure of the disk is small, the 
effect on the SED might be non-negligible. 

2.2. Inner rim 

We assume in the following that the disk is truncated 
on the inside by dust evaporation. The disk may continue 
towards the star in purely gaseous form, but, as long as 
the gas is optically thin to the stellar radiation, we can 
neglect its effects and describe the disk structure as if it 
has a true inner hole. 

The inner disk radius i? r ; m is determined by requiring 
that the disk temperature equals the dust evaporation 
temperature, which we take T ovap = 1500K. The inner 
rim of the disk exposes a vertical surface to direct stel- 
lar radiation and will therefore be hotter than flaring disk 
model at the same radius. For a given inner radius i? r i m 
of the disk, the blackbody temperature of the rim is: 
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(l + ff rim /i? rim ) 1/4 , (11) 



which is derived by equating the emitted blackbody flux 
cT4 m to the received flux from direct stellar radiation 
L»/47ri? 2 im . The factor (1 + H r i m / R^) 1 ^ approximately 
accounts for the effect of self-irradiation, i.e. irradiation 
of the rim by its own emission. The ratio H r i m / i? r ; m is 
an approximate expression for the covering fraction of the 
rim enclosing the cavity R < i? r ; m . We assume i?* <C i? r i m 
and therefore neglect that the star itself obscures a small 
region on the far side of the rim. The vertical height of 
the inner rim, i? r ; m , is given by 
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where /i r ; m is the pressure scale height at the inner rim 
, _ / kTn m R^ im 



4 



and Xrim is a dimensionless constant (similar to %cg de- 
fined above). In appendix A. 3 we show how to compute 
X r i m self-consistently. 



By equating T r \ 
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1500K one can solve for 



-Rrim by writing Eq.(ll) as 
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and iterate until convergence. If the self-irradiation can 
be ignored (H r - im / i? r i m -C 1), then no iteration is needed. 
However, as it turns out, for most HAeBe stars the ratio 
Hrim/Rrim (which is the covering fraction of the inner rim) 
is of the order of 0.1 to 0.25, and therefore the effects of 
self-irradiation cannot be ignored. 

Self- irradiation has the effect of pushing out the inner ra- 
dius, if the temperature of the inner rim is fixed at 1500K. 
This means that the blackbody-emitting surface is larger. 
But since only part of this radiation will escape to infinity 
(the other part being again absorbed by the inner rim it- 
self), energy remains conserved. Once i? r im is determined, 
the vertical height of the inner rim H rim automatically fol- 
lows from Eqs.(12,13), and one has the basic properties of 
the inner rim fixed. 

2.3. Shadowed region 

The inner rim will always be significantly higher than 
the flaring disk height, because it is hotter by a factor of at 
least (2/a) 1 / 4 . In general this will roughly yield a puffing- 
up ratio of H r i m /H cg ~ 2. The disk material immediately 
behind this inner rim will therefore be deprived of direct 
stellar light, since the rim casts a shadow over the disk 
behind it. This region of the disk is cooler than predicted 
by CG97 models, and therefore not flared. Only at much 
larger distance, where the disk comes out of the inner rim 
shadow and is illuminated again by the stellar light, flaring 
may occur. Since the shadow covers every point (R, z) for 
which z/R < H r i m /R r i m , the radius Ra at which the disk 
can again be described by the CG97 model is computed 
by solving the following equation: 

Ra 

-ffrim(-Rrim)-^ = H cg (Rfi) . (15) 

For R < Rd, the most important sources of irradiation 
are the near-infrared emission of the rim itself and the 
fraction of stellar light that is scattered toward the disk 
by the edge of the rim, or by a halo of dust surrounding 
the star+disk system (Natta 1994). In addition to this, 
radiative diffusion can exchange energy between neigbor- 
ing annuli. However, rough estimates of these effects show 
that the radiation emitted by this shadowed region of the 
disk is negligible compared to the emission of the rim and 
of the outer disk, directly heated by the star. 

However, the disk structure immediately behind the in- 
ner rim is of importance for our discussion, because we 
need it in order to estimate the quantity x r ; m in Eq.(12), 
i.e. the ratio of the photosphcric to the pressure scale 
height at i? r i m (see Appendix A. 3). At this location of the 
disk none of the irradiation sources play any significant 
role. Instead, direct radiative diffusion from the inner rim 
will dominate the energy balance there (see Fig. B2). It 
is difficult to determine exactly what the structure of the 
disk is in this diffusive region, but a simple model for this 
will be provided below. 



The inner rim of the diffusive region receives direct 
starlight. This radiation is absorbed and re-emitted by 
the dust grains. Some of this radiation will once again 
be absorbed and re-emitted. In the end, most radiation 
will, after one or more absorption/re-emission events es- 
cape from the hot inner rim surface. But a small fraction of 
the radiation remains trapped within the disk and diffuses 
towards larger radii. The diffusive flux is proportional to 
the gradient of the temperature of the disk's interior. Since 
it is expected that the temperature has a net decrease out- 
wards, there will be an outward pointing diffusive flux. As 
a package of energy has moved one surface height H to- 
wards larger radii, it has had an equal chance of having 
escaped from the disk's surface. This is true for each sur- 
face height it travels outwards. The differential equation 
for the disk's temperature would be approximately: 

d(RT^) _ i?T 4 

dR ~ H ' { ' 

where T is a rough estimate of the temperature of the disk 
interior at radius R. In this equation we have included the 
effect of geometric dilution due to the cylindrical geome- 
try, hence the factor R inside the brackets of the differ- 
entiation. A similar factor H is not included, since H is 
expected to decrease outwards, which would lead to an 
unphysical geometrical compression of the radiation field. 
In this diffusive region the disk could be vertically non- 
isotermal, but we ignore this effect in this simple treat- 
ment. Also we ignore the fact that the definition of H, 
as given in Section 2.1 is not very meaningful in this non- 
irradiated region. Now, since H oc VTR 3 (see e.g. Eq.5, 
and assume for convenience that \ =const) the temper- 
ature T can be eliminated from the differential equation, 
leaving an ordinary differential equation to be solved for 
H: 

rflg(ff 8 ir n ) _ _ 1_ 
dR ~~ H ' 
This equation can be solved numerically. For H <C R, it 
reduces to 



(17) 



d(H/R) 
dR 
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(18) 



which shows that at the inner edge H/R approximately 
goes asH/R= H(R lim )/R lim - \g(R/R lim )/8. In this ap- 
proximation the disk has lost all its energy at R — i? r i m = 
8H(R I i m ), and has collapsed to zero height. This would 
be of course unphysical, and before this happens, other 
heating mechanisms will have taken over, such as the ring 
emission mentioned above, light scattering off a spherical 
halo surrounding the disk+star system, interstellar radia- 
tion or viscous dissipation. 

Eq.(17), but with an opposite sign, describes the diffu- 
sion of energy inwards from the still flaring part of the disk 
into the shadowed region. 

2.4. Self-irradiation corrections to flaring part 

For R> Rfi the flaring disk is in full view of the central 
star, and the irradiation is clearly dominated by the stellar 
flux. However, for HAeBe stars the luminosity of the in- 
ner rim can be of order 20% of the stellar flux. Moreover, 
the emission of the inner rim is not isotropic. In close to 
polar directions, the projected surface of the rim is small 
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in the flux correspondingly low. For a point on the sur- 
face of the flaring disk on the other hand, the projected 
surface of the rim is close to 4_R r i m if r ; m (see Appendix B). 
Therefore, the relative flux contribution of the rim at the 
disk surface will even be larger than the angular average 
of 20%. For an accurate determination of the disk struc- 
ture in the flared region one should include this secondary 
irradiation. We consider the inner rim as a cylinder. The 
irradiative flux is then (see appendix B): 



F irr =a(R)F* + -a{R)aT I 
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rim 



R, 



it \ R 

cos 9 [6y/l — 8 + arcsin S] . 



(19) 



This irradiative flux F" r is then equated to the disk's ree- 
mission flux F cm ; t (Eq.3), thus yielding the corrected value 
of Ti. 

For the surface layer temperature one also has to add 
the two sources of irradiation. But in order to compute 
the e s factor in Eq.(7), one now has to replace the Kp(T*) 
in the definition of e s (Section 2.1.2) with: 



Kirr 



(20) 



3. COMPUTING THE SED 

Once the disk structure is known, one can compute the 
emerging spectrum. In principle one requires a 2-D (axi- 
ally symmetric) ray-tracing program to compute the SED 
at various inclination angles. This is because the disk 
structure, in particular for HAeBe stars, is far from be- 
ing flat, and self-occultation is a real possiblity (Chiang & 
Goldreich 1999). Yet, in this section we shall use a simpli- 
fied treatment of inclination, therefore limiting ourselves 
to inclinations that are not too edge-on. 

3.1. Spectrum of flaring disk 

The determination of the spectral components from the 
flaring part of the disk is relatively straightforward, and 
has been described fully by CG97. We confine ourselves 
here to a quick reference of the equations. The emission 
from the interior of the disk is: 



Fl(d,i) = 2?rcos 



[1 — exp(— H{R)k„/ cosi)} 
Rn " (21) 

x B v (Ti(R))^dR, 

where d is the distance from the observer to the source, 
and i the inclination angle (i — means face-on). This 
equation remains valid at wavelengths where the disk inte- 
rior becomes optically thin. The emission from the surface 
layer is 

/■flout 

F°(d,i) = 2n / [1 + exp(-E(i?)/c 1/ /cosi)] 

jRn ' ' (22) 

x B v {T s {R))/^K v ^dR, 

where AS is given by Eq.(10). Eq.(22) includes the effect 
of also seeing the surface layer on the other side of the disk 
at wavelengths where the disk is optically thin. 



3.2. Spectrum of inner rim 

The inner rim is assumed to be a cylinder emitting on its 
inner side as a blackbody with temperature T r ; m . To com- 
pute the emission from this rim we have to include inclina- 
tion effects. Clearly, if i = (facc-on), then the emission 
from the rim is zero (although in reality the rim may be 
more rounded-off, and therefore also visible at face-on in- 
clination). For i — it/2 (edge-on) the emission is also zero, 
because the rim is self-occulting. The projected surface, 
as seen by an observer at inclination i is given by Eq.(B2). 
For not-too-large inclination angles we are in the regime 
6 <C 1 and one can then write: 

Fr(d,i) = 4 ^f rim sin^(r rim ) • (23) 

The total luminosity of the rim equals L r ; m = ri r ; m L + , 
where f2 r i m = i? r im/-Rrim is the covering factor of the rim 
(i.e. the fraction of the sky covered by the rim as seen from 
the perspective of the star). For large H^ m /R rim this will 
therefore give a strong contribution to the SED. 

4. RESULTS 

In this section we show the structure of the computed 
models and the SEDs produced by models with different 
parameters. 

4.1. Example model 

To illustrate the general features of our model we set up 
an example model consisting of a star with T e s = 9500K, 
L* = 47 Lq and M* = 2.4 M Q (we take the stellar spec- 
trum to be a blackbody for simplicity), and a disk with 
outer radius i? ut = 400 AU and inner radius at the dust 
condensation point (T rim = 1500K), which for this setup 
is at i? r i m = 0.47AU. The disk's surface density goes as 
E(i?) = 2 x 10 3 (R/AV) , with /3 = -1.5 and we assume 
a gas to dust ratio of 100. The disk is seen at an inclina- 
tion angle i = 45°. We take a simple model for our dust 
opacity: a mixture of amorphous carbon and astronomical 
silicate in a ratio 0.05. The amorphous carbon opacity is 
computed according to the simple recipe given by Ivezic 
et al (1997), while the silicate opacity is that of Draine & 
Lee (1984). We ignore scattering. The SED depends on 
dust properties in various ways, and variations of the opac- 
ity from object to object and between grains in the disk 
surface layer and midplane arc known to occur (e.g. Bouw- 
man et al. 2000). However, this aspect of the problem is 
outside the scope of this paper, and is not relevant to the 
conclusions we reach. 

Fig. B3 shows the structure of the disk for these stel- 
lar parameters. The inner edge is at R = 0.52AU, beyond 
which the vertical height of the inner edge quickly declines 
(although when plotted linearly in R this decline is not so 
rapid as it seems in this figure). For R > 6.6AU the disk 
has the usual flaring shape , but with self-irradiation ef- 
fects added. The region between 0.52AU < R < 6.6AU 
is the shadowed region. Between 0.52AU < R < 1AU 
the disk is dominated by radiative diffusion from the inner 
edge outwards, as described in Section 2.3. In the region 
1AU < R < 1.5AU radiative diffusion from outside in- 
wards dominates. For radii 1.5AU < R < 6.6AU the disk 
is still not in sight of the star itself, but it receives flux from 
the inner rim, which is sufficiently strong to keep up the 
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disk. Inwards radiative diffusion prevents a sudden jump 
to the fully flared region starting at R > 6.6AU, although 
the surface temperatures cannot be prevented from jump- 
ing. The inward radiative diffusion is estimated using an 
equation similar to Eq. (16), but with opposite sign. The 
details of the structure between 1AU < R < 6.6AU are 
not very important, since this region hardly contributes to 
the SED. For this reason we are content with the rather 
crude description presented here. 

Fig. B4 shows the resulting SED for the example model. 
We show for reference in panel (A) the SED predicted by 
a CG97 model, where the effects of the inner radius (emis- 
sion and shadow) is ignored. Panel (B) shows the effect 
of including the emission from the inner rim, assumed to 
have the height of the flaring disk at that radius. In this 
case, the SED in the mid to far infrared remains unmod- 
ified, but a new component appears in the near-infrared. 
Finally, panel (C) shows the results of the self-consistent 
calculation described in the previous sections. One can see 
that the NIR emission from the puffed-up rim increases, 
but the inner rim starts to cast a shadow over the disk, 
thereby suppressing some of the mid- and far infrared emis- 
sion. Panel (D) shows a blow-up of the 1-100 /xm region, 
where the different effects of the inner rim can be clearly 
seen. 

In the following subsections we will show how the model- 
predicted SEDs change as function of some of the model 
parameters. We summarize some of the characteristics of 
the models and of the SED properties in Table 1. 

4.2. Effect of inclination angle 

In Fig. B5 we show how the SED of our example model 
changes with the inclination with respect to the line of 
sight (Panel A). We consider here only values of i such 
that the line of sight does not intercept the outer disk. 

One sees that for nearly face-on inclinations the NIR 
excess is less prominent and the far infrared is strong, 
while for nearly edge-on inclinations the opposite is true. 
By varying i from 50° to 20°, the ratio Fl ^ ,R decreases 
from 0.18 to 0.08, while the infrared excess at wavelengths 
longer than 7 /im increases from 0.32 to 0.35. This is be- 
cause the disk radiates predominantly towards the poles, 
while the inner rim radiates predominantly towards the 
equator. The emission features are optically thin features 
and are therefore less affected by the inclination. Only 
when the inclination would exceed the opening angle of 
the flared disk, these features would go from emission into 
absorption because of the obscuration effect of the disk 
itself. However, in order to model this properly we would 
need a multidimensional ray-tracing computation, which 
is beyond the scope of this paper. 

4.3. Effect of varying inner radius 

The inner rim in our models is located at the dust evap- 
oration radius, and is therefore fixed once the stellar pa- 
rameters are defined. Nevertheless, it is instructive to see 
the effects of varying inner radius on the SED, because it 
helps in clarifying the properties of our models, and we 
plot in Fig. B5 the SED of our standard model for varying 
i? r i m (Panel B). The effect is strong in the near-infrared, 
where both the excess luminosity and its shape change 
drastically. For larger i? r ; m , the near-infrared "bump" is 



enhanced, with a peak that moves toward longer wave- 
lengths, but the flux around 1/im decreases as the temper- 
ature of the inner rim decreases. 

There are also, barely visible in the figure, some changes 
at intermediate wavelengths, due to the reduction of shad- 
owing which increases the contribution of the disk surface 
layer also around 50/xm. 

4.4. Effect of shadow on the 10 micron feature 

If the height of the puffed-up inner rim is computed 
according to the equation for vertical pressure balance 
(Eq. (5) with 7] replace by T r i m ), then the shadow that 
the rim projects over the disk will usually extend up to 
a radius of roughly 10 to 20 times the inner radius. This 
ratio turns out to be fairly robust in our model, which 
can be attributed to the fact that the ratio of T r ; m /T cg 
is fairly constant for all model parameters. At 10 to 20 
times i? r ; m the surface layer temperature is about 4.5 to 
6.3 times lower than T r j m (which is 1500K by definition in 
our model). This means that the shadow starts to affect 
the regions of the disk's surface layer that emit the 10 mi- 
cron feature, which lies at 5 times longer wavelengths than 
the NIR bump. This effect is already visible when com- 
paring our standard model to the result of a CG97 model 
(Fig. B4C and D). 

If, for some reason, we have underestimated the height 
of the inner rim (which is the most difficult quantity to 
determine in our models), this might have repercussions 
on the predicted 10 micron flux. To investigate this effect, 
we have artificially enhanced the height of the puffed-up 
inner rim (i.e. we have relaxed the assumption that the 
disk is here in vertical pressure balance), and computed 
the SED. The results are shown in Fig. B5C, where we 
have increased \nm to values as large as 8. For increasing 
inner rim height, the 2 micron bump is of course enhanced, 
but the major effect on the SED is the change at longer 
wavelengths, due to the fact that the shadow extends to 
larger radii, eventually to most of the disk. The 10 micron 
feature is the first part of the SED to be affected by the 
increase of Xrim, and has entirely disappeared for Xrim=8. 
The 20 micron silicate feature is less affected, because it is 
emitted at lower temperature, and hence at larger radius. 
The total IR flux is conserved: the increase in NIR flux 
goes at the cost of the mid-IR and far-IR flux (see Table 
1, models H1..H4). 

The exact value of Xrim depends on a number of prop- 
erties that are difficult to constrain, such as the shape of 
the disk behind the rim and the dust properties. However, 
its exact value is not crucial for the discussion of the basic 
idea behind our model, which revolves around two issues: 
the presence of a puffed-up inner rim at roughly the dust 
evaporation radius, naturally producing the NIR bump, 
and a proper consideration of energy conservation, which 
implies that we must take into account the shadow the rim 
projects onto the disk. The exact shape of the rim is, in a 
way, a secondary aspect. 

4.5. Effect of Y,-slope 

The surface density T,(R) is an unknown quantity, and 
serves as an input to the model. The results of the models 
turn out not to be very dependent on the surface density. 
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The main effect it can have is changing the balance be- 
tween near-IR and far-IR radiation, and the slope of the 
SED at far-IR/submm wavelengths. To demonstrate this, 
we show models with E = E (i?/-Ro) /3 for different values 
of [3. For each of these models we tune the Eo such that 
the total covering fraction of the disk Q = H(R out )/R out 
is equal to ft — 0.4, see table 1 for the resulting values of 
E . The rest of the parameters are the same as for the 
example model. The results are shown in Fig. B5D. 

4.6. Behavior of the model for varying stellar type 

The behaviour of the near-infrared excess as function 
of the stellar type of the central star is shown in Fig. B6, 
which plots the SEDs of four stars along the zero-age-main- 
sequence ranging from spectral type B2 to M2 (see Table 
1). The NIR bump is clearly visible in stars of spectral 
type A and earlier, but tends to disappear in later stars. 
This is due to two effects. First of all, for cooler cen- 
tral stars the stellar flux peaks at longer wavelength, and 
therefore swamps the near-infrared emission of the disk. 
In addition to this, the ratio of the covering factor of the 
puffed-up inner rim (in vertical hydrostatic equilibrium) 
to the covering fraction of the flaring part of the disk de- 
creases for later spectral types. These factors determine 
the relative importance of the various components of the 
SED, as can be seen from Tabic 1. As we consider stars 
of later spectral types, the near-infrared excess luminos- 
ity decreases by a factor of two. The various component 
that contribute to the observed SED of a typical TTS are 
shown in Fig. B7. 

Our disk models provide a natural explanation for the 
fact that the near-infrared "bump" is observed in Herbig 
AeBe stars, while it is not noticeable in stars of later spec- 
tral type, and in particular in TTS. We will come back to 
this point in the following section. 

5. DISCUSSION 

5.1. Disk physics 

It is currently not clear whether a sharp inner rim can 
be stable in the way described here. The thermal pres- 
sure which is responsible for the rim height will act not 
only in the vertical but also radially toward the star. We 
can think of two possibilities to stabilise the rim. The 
space between the star and the disk could be filled with 
optically thin gas, and the gas pressure could stabilize the 
rim. A second possibility would be centrifugal force. If 
the inner rim rotates at slightly super-keplerian speeds, it 
could counter the gas pressure. This super-keplerian ro- 
tation might be automatically created if material flowing 
inwards because of an initial pressure gradient preserves it 
angular momentum. A dynamical study of the problems 
would be very interesting, but is beyond the scope of the 
current study. 

Radiation pressure on dust grains may also affect the 
structure of the inner rim. In the upper regions, the 
gas drag may not be strong enough to keep small dust 
grains from being blown away by radiation pressure, an 
effect which would locally decrease the dust-to-gas ratio 
and therefore opacity and the rim height. On the other 
hand, radiation pressure could also lead to gas-dust sepa- 
ration at the inner rim if some low level accretion is still 
going on in the disk. The star would then accrete only 



the gas while the dust would stay in the rim and accumu- 
late with time. This might cause an increase of the rim 
height, compensating for the removal of grains by radia- 
tion pressure. Radiation pressure on dust grains can also 
have effects in the flaring part of the disk, by pressing the 
grains deeper into the disk and reducing the surface height 
there. All these processes are poorly studied and deserve 
more attention in future work. 

Another point of concern is the question whether the 
flaring part of the disk is dynamically stable. If one takes 
the disk equations of CG97, and those described in the 
present paper, at face value, one can show that the disk 
is unstable to self-shadowing effects caused by a perturba- 
tion in the vertical scale height (Dullemond 2000, Chiang 
2000). In fact, the shadowed region outside the inner rim 
is a an example of this shadowing effect. This instability 
acts on a Keplerian time scale and could cause the entire 
flaring disk to collapse into the shadow of the inner rim. A 
recent study however shows that multi-dimensional radia- 
tive transfer effects can stabilize the disk for certain disk 
parameters (Dullemond submitted). And the inclusion of 
self-consistent vertical structure of the disk models may 
enhance this stabilizing effect. 

5.2. Comparison with the observations 

Disk models where the irradiation from the star domi- 
nates the heating provide a good fit to most properties of 
pre-main-sequence stars over a large range of masses (see, 
for example, D'Alessio et al. 1998, 1999, CG97, Chiang 
et al. 2001). However, there are few points that have re- 
mained unclear, and we think that the better treatment of 
the inner disk with respect to previous models provides an 
answer to one of them, namely the long-standing puzzle of 
the near-IR bump in the SEDs of HAcBe stars. 

In Fig. B8 we show a fit of our model to the SED of 
AB Aurigae, one of the best studied HAe stars. The basic 
features of the SED arc well reproduced by the model: the 
3/im bump, the 10 micron feature, the far IR plateau and 
the sharp decline towards mm wavelengths. Our model 
has only a few free parameters (z, E(i?), R out and the dust 
opacity table), so the fact that the data can be fit in a rea- 
sonable way gives support to our model. The inclination 
angle of i = 65° needed for this fit is larger than what most 
observations suggest (e.g. Natta et al. 2001 and references 
therein). However, the assumption of a perfectly vertical 
rim which is responsible for this value may be too simple. 
A more rounded surface of the inner rim will allow more 
of the rim radiation to be emitted into polar directions. 

Nevertheless, the model-predicted values of F p IR for A 
stars (hovering around 0.17) agree well with those ob- 
served in four HAe stars by Natta et al. (2001), which 
range between 0.12 and 0.25, respectively, and by Meeus 
et al. (2001) in a sample of 14 HAe stars, where 
varies between 0.06 and 0.3 for the group I sources. 

Observations show that the "bump" is a more or less 
constant feature from star to star (Hillenbrand et al. 
1992), in contrast to the long- wavelength emission, which 
is affected by a variety of disk and dust parameters (Meeus 
et al. 2001). This robustness is naturally reproduced in 
our model, since to first order the properties of the NIR 
bump are determined only by the dust sublimation tem- 
perature and the covering fraction of the rim, the latter 
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varying only weakly as a function of stellar type. 

In spite of this, the "bump" disappears in stars of spec- 
tral type F or lower (see Fig. B6). As discussed in §4.6 
and shown in Table 1, this can be explained as being due 
to two concurring effects. The first is the lower covering 
angle of the rim in stars of later spectral type, where dust 
sublimation occurs closer to the star. The second is the 
fact that the stellar radiation peaks at longer wavelengths, 
swamping the rim emission. Thus, even if the shape of the 
"bump" is the same, it becomes very difficult to see it and 
to separate the rim contribution from the stellar emission. 
It is, however, possible that one is seeing the rim emis- 
sion in some TTS with near-infrared colors that can only 
be understood if the disk has an inner radius of few stel- 
lar radii, coinciding with the dust sublimation radius for 
systems of typical TTS luminosity (Beckwith et al. 1990; 
Kenyon et al. 1996). 

Another interesting aspect of these models is the effect 
of the shadow that the rim projects onto the disk. The 10 
fim silicate emission in our example model has -^-=3.5, 
in reasonable agreement with what observed in HAe stars 
by Natta et al. (2001). These numbers, however, are very 
dependent on the properties of the grains in the disk sur- 
face layer. A weak or absent 10 fim feature, for example, 
can be attributed to a dominance of large silicate grains in 
the disk surface layer (Meeus et al. 2001), to a very high 
ratio of the visual opacity to the opacity in the feature 
(Natta et al. 2001), or to an almost perfect cancellation 
of absorption and emission in disks seen edge-on (Chiang 
& Goldreich 1999). We propose here still another possi- 
bility, namely that the inner rim in some object is higher 
than normal, because, for example, the rim is not in hy- 
drostatic equilibrium, is thicker in the radial direction, or 
has a higher Xrim, the ratio of the photosphcric to the pres- 
sure scale height. As shown in Fig. B5C and Table 1, a 
factor 1.7 increase of Xrim would be sufficient to suppress 
the silicate feature entirely. 

One point we want to stress here is that a rim structure 
such as we have described can only exist if the gaseous disk 
is optically thin to the stellar radiation and the dusty disk 
is optically thick. If an optically thick gaseous disk extends 
toward the star closer than the dust evaporation radius, 
then there is no puffed-up rim and the standard models, 
which describe the emission of an optically thick disk sur- 
face, apply. The formation of a puffed-up rim seems pos- 
sible only if the accretion rate through the disk is very 
low, as briefly discussed by Natta et al. (2001). There are 
indications that this is indeed the case in the majority of 
HAe stars. Tambovtseva et al. (2001) find accretion rates 
below 10~ 8 Mq yr _1 from their study of Balmer emission 
lines in UX Ori (see also Ghandour et al. 1994). There 
are no similar studies for the hotter HBe stars, and it is 
clear that in many TTS accretion rates are high enough to 
ensure that the inner gaseous disk will be optically thick 
(Gullbring et al. 1998). This is an additional explana- 
tion why the near-infrared "bump" is seen only in a minor 
fraction of TTS, if at all. 

The idea that the emission in the near-infrared is domi- 
nated by the emission of the inner rim finds support in the 
recent intcrfcrometric observations of Millan-Gabet et al. 
(2000). They find that in J and K the emission of AB Aur 
is best fitted by the emission of a thin shell of radius ~ 0.3 



AU, almost spherical on the plane of the sky. For the lu- 
minosity of AB Aur, the inner rim is located at 0.5 AU 
for dust evaporation temperature around 1500 K. For the 
low inclination derived from millimeter and optical data 
(i<45°; see Natta et al. 2001 and references therein), the 
rim projected on the plane of the sky is almost spherical. 
We expect that other HAe stars will be resolved by current 
and future interferometers, so that the main feature of our 
models, namely the fact that the near-infrared emission 
comes from a narrow rim located at the dust evaporation 
radius, will be directly checked. 
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APPENDIX 
FINE-TUNING FUNCTIONS 

The tp parameters 

When the interior of the disk is not optically thick anymore to its own radiation and/or the radiation from the surface 
layer, then ip will deviate from unity. For an optically thick interior, the total emitted flux (one a single side) is simply 
-Femit = <jT*. For non optically thick cases this becomes 



/>oo 

Fomit - / 7rB v (Ti) [1 - cxp(-£^)]di/ 
Jo 



(Al) 
= ipi oT^ . 

with 

_ S"B V (TQ [1 - expj-EK^du 
Jo B v {Ti)dv 

The limits for high and low optical depth are 

f 1 for S« P (T i )»l . . 

S Kp (Ti) for £ Kp (Ti) < 1 [Ad) 

The absorbed flux from the surface layer equals the flux emitted downwards by the surface layer multiplied by the 
absorption fraction [1 — exp(— £«„)], and integrated over frequency: 



poo 

F abs = / 2ttA£ k v B v {T s ) [1 - exp(-£K„)]di/ 
Jo 



(A4) 

where F s is the total flux from the surface layer downwards, and 

B v (T s )k v [1 - exp(-EK v )]di> 
4> s = — r oo p (rr : , (A5) 

The limits for high and low optical depth are 

/ _/ 1 for £ KP (T S )»1 , AR x 

^ s -\ E4(T S )/ KP (T 3 ) for S KP (T S )«1 ^ 



where kq is the Planck square mean opacity: 

,^ (T) , (A7) 



2 = S~B v {T)nld V 
Q[ >~ $™B v {T)dv 



To a certain degree of accuracy the approximation Kq(T s )/kp(T s ) ~ «p(T s ) can be used. 

When the above limiting expressions for (Eq.A6) and ip\ (Eq.A3) are inserted into Eq.(4), then one reproduces 
Eqs.(12a,12b,12c) of CG97. 

Since these f/'i and ip s depend on 7] itself (the determination of which requires the values of ipi and tp s ), the solution of 
Eq.(4) requires an iterative procedure. 

The Xcg parameter 

For a given grazing angle a, the ratio x cg = H/h is independent of h. However, a itself depends on H through Eq.(l). 
We define H to be the height at which 



S f p(z)dz = 1 (A8) 
" Jh 



Since the vertical density distribution is assumed to be a Gaussian, this amounts to solving the following equation 

where, for a given value of h, the a depends on % cg through Eq.(l). Eq.(A9) can be quickly solved for \cg using any kind 
of root-finding algorithm. 
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The Xrim parameter 

In this appendix we wish to determine the dimensionless parameter Xrim = ^rim/^rim, which determines the vertical 
height of the inner rim. We define the surface height to be the height to which the optical depth of the rim on a radially 
outward directed ray is greater than 1. 

/>oo 

T r {z )= p(R,z)n P {T,)dR=l, (A10) 

where z — zqR and z is the value of z at R = i? r ; m . We neglected a geometric factor a/ 1 + Zq / R^ im (arising from 
cylindrical coordinates) because it is very close to unity. 

In order to determine r r (zo) we need to know the two-dimensional structure of the puffed- up inner rim, which we 
discussed in section 2.3. We have estimated the radial behavior of H{R)/R to be linear with a slope of —1/8. A radial 
ray through the upper layers of the inner edge will therefore have a roughly 8 times higher optical depth than a vertical 
ray between z = z and z = oo. In order to estimate zo we therefore estimate r r (z ) to be simply 8 times t v (z ): 

/oo 
p(i? rim , z)K P (T*)dz' 

(All) 

e x dx , 



8£k p (T*) 



which can be easily computed following the procedure in appendix A. 2. 

EMISSION FROM A CYLINDER 
Far field limit 

Consider the inner rim to be a cylinder with radius i? r i m and vertical height H r i m , emitting toward the inside as a 
blackbody of temperature T r ; m . The unocculted surface as seen by an observer at inclination i (measured from the pole) 
can be computed as follows. Define the quantity 6: 

S=^ t <mi. (Bl) 

5 is defined so that for 6 = 1 the inclination is just small enough that we are able to see the central star (which is assumed 
to be a point-source). For 6 < 1 we have a more face-on view, while 5 > 1 means a more edge-on view of the system. The 
observed flux at a distance d from this cylinder in the far-field limit is for S < 1 

2 



F v = 2B v (T Tim ) (^^j cosi SVl-S 2 + &TCsmS , (B2) 



and for 5 > 1 



Rr 



F v =wB v (T lim )\-^-j cos*. (B3) 
Near field limit 

In the far-field limit, the rim becomes invisible at i = n/2 which corresponds to z/R = 0. In order to compute the 
contribution of the inner rim to the irradiation of the flared disk, we need to correct the equation for the fact that the 
rim becomes invisible already at z = H rim . We do this by introducing a modified inclination 9 defined by 

tanfl= ^~g" m (B4) 

where we took the reference point to be the nearest corner of the cylinder instead of the star. We can then write 5 as 

= R/R vim - 1 

z/H rim - 1 1 ' 

The flux is then further given by Eqs.(B2,B3). Comparing the resulting flux with exact numerical evaluation of the 
projected surface has shown that this approximation gives irradiation fluxes accurate to within 20% for inclinations > 7r/4 
and R > 2i? rim . 
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Figures 




Fig. Bl. — A sketch of the geometry of the model. 
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Fig. B2. — A sketch of the two dimensional radiative diffusion process through the disk inner edge. The solid arrows indicate direct radiation 
from the star. The dashed arrows represent the diffusion through the disk. 
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Fig. B3. — The structure of example model E3. Top panel: disk pressure scale height (solid line) and surface height (dashed line), both 
divided by radius. Bottom panel: midplanc temperature (solid line) and temperature of the surface layer (dashed line). 
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Fig. B4. — SED of the example model (see section 4.1), computed according to the unmodified CG97 model (panel A, model code El), 
the CG97 model with inclusion of the emission from the inner rim (panel B, model code E2), and with a self-consistent puffed-up inner rim 
(panel C, model code E3). In panel D the three SEDs are compared. In the unmodified CG97 model, the SED (solid line) is the sum of the 
emission of the star (dash-dotted line), the disk interior (dashed line) and the disk surface (dotted line). When the emission from the inner 
edge is added, without modifying the vertical height of the inner rim, then a blackbody component at NIR wavelengths is added to the SED 
(long-dashed line in Panel B). When the height of the inner rim is consistently computed, the SED is also the sum of four components (same 
symbols as in Panel B), but the NIR emission increases, and a shadow is cast over the disk, thus suppressing some of the disk emission. 
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Fig. B5. — The observed SED of our example model, for varying input parameters. For all figures the inclination angle is i = 45 degrees, 
except for the one that shows the effect of varying inclination. 



17 




Fig. B6. — Model-predicted SED for ZAMS stars of different spectral type. The stellar parameters are given in Table 1. 
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Model applied to T-Tauri disk : 




0.1 1.0 10.0 100.0 1000.0 

A (fim) 



Fig. B7.— SED of a T Tauri star of T eS = 3800 K, L/Lq = 0.76 and M/M e = 0.55. The various line styles represent the different 
components contributing to the SED in the same way as in Fig.(B4). 
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Fig. B8.— Best fit to AB Aur. Parameters are M« = 2.4M Q , L, = 47L Q , T* = 9520K, E = 10 4 * (i?/AU)~ 2 , R out = 400AU, i = 65°, 
D = 144pc. The plot shows the model fit (solid line), together with photomctcry points (squares) and the ISO SWS and LWS spectra between 
2.3 and 200 /im (data from Meeus et al. 2001). 
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Tables 



Table Bl 



Properties of different models discussed in this paper. 



Scries 


Fig 


Model 


T c ff 


F, 


So 







Rn 


Hrim 


Xrim 


-Fnir 


Fir 




Fio 


F& 


gem- 2 


AU 


AU 


flrim 


F, 


F, 


Fir 


F7.7 


Example 


B4 


El 


9500 


47 


2000 


-1.5 


0.47 


0.0 


0.11 


5.3 


0.03 


0.35 


0.08 


4.7 






E2 


9500 


47 


2000 


-1.5 


0.50 


0.0 


0.10 


2.9 


0.12 


0.50 


0.24 


3.6 






E3 


9500 


47 


2000 


-1.5 


0.52 


6.6 


0.18 


5.0 


0.17 


0.50 


0.33 


3.5 


Inclination 


B5A 


11(50°) 


9500 


47 


2000 


-1.5 


0.52 


6.0 


0.18 


5.0 


0.18 


0.50 


0.36 


3.2 






12(35°) 


9500 


47 


2000 


-1.5 


0.52 


6.0 


0.18 


5.0 


0.13 


0.47 


0.28 


3.9 






13(20°) 


9500 


47 


2000 


-1.5 


0.52 


6.0 


0.18 


5.0 


0.08 


0.43 


0.19 


5.2 


^rim 


B5B 


Rl 


9500 


47 


2000 


-1.5 


0.20 


2.2 


0.15 


5.3 


0.10 


0.48 


0.21 


6.6 






R2 


9500 


47 


2000 


-1.5 


0.40 


3.5 


0.17 


5.1 


0.15 


0.50 


0.30 


4.4 






R3 


9500 


47 


2000 


-1.5 


0.80 


7.9 


0.20 


4.9 


0.18 


0.50 


0.36 


2.3 




B5C 


HI 


9500 


47 


2000 


-1.5 


0.52 


6.0 


0.18 


5.0 


0.16 


0.49 


0.34 


3.4 






H2 


9500 


47 


2000 


-1.5 


0.52 


17.0 


0.22 


6.0 


0.20 


0.48 


0.43 


1.9 






H3 


9500 


47 


2000 


-1.5 


0.53 


35.0 


0.26 


7.0 


0.25 


0.49 


0.50 


1.2 






H4 


9500 


47 


2000 


-1.5 


0.54 


70.0 


0.30 


8.0 


0.29 


0.48 


0.60 


0.8 


S power 


B5D 


SI 


9500 


47 


5 


-0.5 


0.50 


2.9 


0.12 


3.5 


0.11 


0.48 


0.23 


4.8 






S2 


9500 


47 


100 


-1.0 


0.52 


4.2 


0.16 


4.3 


0.14 


0.48 


0.29 


3.9 






S3 


9500 


47 


2000 


-1.5 


0.52 


6.0 


0.18 


5.0 


0.16 


0.49 


0.34 


3.4 






S4 


9500 


47 


4T0 4 


-2.0 


0.52 


6.0 


0.20 


5.6 


0.19 


0.51 


0.37 


3.2 


ZAMS 


B6 


Z1(B2) 


22000 


1000 


2000 


-1.5 


5.88 


68.0 


0.27 


4.5 


0.31 


0.51 


0.60 


2.6 






Z2(A0) 


8900 


50 


2000 


-1.5 


0.53 


6.1 


0.18 


5.0 


0.16 


0.48 


0.34 


3.4 






Z3(G2) 


5900 


1 


2000 


-1.5 


0.08 


0.73 


0.12 


5.4 


0.10 


0.48 


0.21 


3.3 






Z4(M2) 


3500 


0.05 


2000 


-1.5 


0.015 


0.22 


0.09 


5.7 


0.07 


0.48 


0.15 


2.2 


T Tauri 


B7 


Tl 


3800 


0.76 


2000 


-1.5 


0.064 


0.64 


0.14 


5.4 


0.13 


0.60 


0.21 


3.0 


AB Aur 


B8 


Al 


9520 


48.3 


10 4 


-2.0 


0.52 


8.1 


0.19 


5.4 


0.22 


0.45 


0.49 


2.5 



The first three columns are for model identification, the next seven list model parameters and disk properties, and the last 
four columns list computed observables evaluated from the star-subtracted SED, where -Fnir is the flux in interval 1.25-7 
/zm, Fir the flux at A > 1.25 /im, and F 10 and F7.7 the fluxes at 10/xm and 7.7/im, respectively. F* is the integrated stellar 
flux. 



